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ОПРЕДЕЛЕНИЕ РЕЗОНАНСНЫХ ДИАГРАММ ЛОПАСТЕЙ 
МАЛОГАБАРИТНЫХ ВЕТРОЭНЕРГЕТИЧЕСКИХ УСТАНОВОК 


А.С. ГУРИНОВ, В.В. ДУДНИК 
(Донской государственный технический университет) 


Представлена методика определения резонансных характеристик лопасти малогабаритной ветроэнергетиче- 
ской установки. Для ветротурбины диаметром 4,5 м были построены расчетные формы первых трех тонов 
колебаний лопастей и построена резонансная диаграмма. 

Ключевые слова: ветроэнергетическая установка, резонансная диаграмма, формы колебаний. 


Введение. Источники возобновляемой энергии — один из современных способов получения элек- 
троэнергии с использованием природных явлений. Ветроэнергетические установки среди прочих 
нетрадиционных источников получили самое широкое распространение. Ветроэнергетика вызы- 
вает повышенный интерес по целому ряду причин. Главной из них является независимость ветро- 
генератора от общей системы электроснабжения, отсутствие выбросов. При правильном выборе 
ветрогенератора его использование приносит существенную экономию средств. Это особенно ак- 
туально при долгосрочной эксплуатации. 

Наилучшие характеристики показывают ветроэнергоустановки (ВЭУ), обладающие систе- 

мами автоматического управления. В случае использования системы автоматического управления 
с нагрузкой необходимо добиваться изменения частоты вращения при частотах, близких к резо- 
нансным. Только в этой ситуации можно добиться длительной работы ВЭУ. Наиболее опасны ре- 
зонансные колебания лопастей в направлении перпендикулярном плоскости вращения. Частоты 
таких колебаний могут быть определены аналитически. 
Методика построения резонансных диаграмм. Ветротурбина, создавая крутящий момент на 
валу генератора, создает и переменные нагрузки, обусловленные как наличием косой обдувки, 
так и дисбалансом самого ветроколеса. Амплитуда переменных нагрузок зависит от конструктив- 
ных особенностей лопасти и режимов работы ветроэнергетической установки. Кроме того, ампли- 
туда колебаний лопастей зависит от того, насколько близки частоты внешних сил, возбуждающих 
колебания к собственным частотам лопастей. При совпадении частот происходит резкое увеличе- 
ние амплитуды, обусловленное резонансом. Для определения резонансов ветротурбины необхо- 
димо построить резонансную диаграмму колебаний в плоскости поперечной плоскости вращения. 

В настоящее время существует ряд методов расчета резонансных характеристик лопастей. 
Однако для малых ВЭУ, самой нижней ценовой категории, необходим наиболее простой метод, 
расчет которого представлен ниже. 

В общем случае уравнение свободных колебаний лопасти в этой плоскости имеет вид [2] 


(ЕЛ „у — (М ву, +т,у, =0, (1) 

где Е, — модуль упругости материала лонжерона лопасти; Г, — момент инерции сечения лонже- 

рона лопасти относительно горизонтальной ОСИ; у = перемещение лопасти в вертикальном 
направлении; №, — центробежная сила; т, — погонная масса лопасти. 


Штрихами обозначено дифференцирование по радиусу, а точками — дифференцирование 
по времени. 

При использовании жесткой заделки лопастей можно принять следующие граничные ус- 
ловия 


Уно) = 0, ре =0, (Е =) в =0, (Е и Уи) в =0. 
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Решение уравнения (1) можно представить в виде произведения двух функций 
У, (",1) =5(1)-6(1), (2) 
где &(7) — функция только радиуса; 65(#) — функция только времени. 
Соответственно выражение распадается на уравнения 


(+ р. 6) =0, (3) 
где р угловая частота собственных изгибных колебаний лопасти. 
[Е78"(^)] - [М5] - рт.) =0. (4) 


Для решения этих уравнений наиболее распространенным является способ Б.Г. Галеркина 
[3], согласно которому решение ищется в виде ряда 


5) =>С6, (5) 


где С, - постоянные коэффициенты; ф,(г) — заранее заданные функции радиуса, удовлетворяю- 
щие граничным условиям; п — количество членов ряда, принимаемых к расчету. 

Функции ф,(г) выбираются так, чтобы каждая из них удовлетворяла всем граничным ус- 

ловиям, а вся совокупность функций ф,(г) сходилась к какому-то значению, а не расходилась. 


Можно воспользоваться функциями Крылова при соответствующем выборе коэффициентов. 
Согласно методу Галеркина ряд (5) подставляется в равенство (4), каждый член умножа- 
ется на одну из функций ф, (г) и интегрируется по радиусу лопасти 


>С Ел, ХС [ибп =0, (6) 


1=00 
где г - номер функции в ряде (5); А- номер функции на которую умножается уравнение (4). 
Для вычисляемых интегралов можно ввести обозначение 


К 

[СЕ „ф’) "фь 4’ =: Ск 

0 

К 

ФУ, (7) 
0 


В 
[т,.ф фа» - Щи 
0 


Таким образом, вместо уравнения (4) можно перейти к системе из п алгебраических урав- 
нений вида (6). Система состоит из бесконечного количества уравнений, но на практике целесо- 
образно ограничиться 4-5 членами и, соответственно, 4-5 уравнениями, которые выглядят сле- 
дующим образом 


Ус (с =8:= роту ) =0. (8) 
20 


Форма колебаний в этом случае определяется с точностью до неопределенного множите- 
ля с›. Этот множитель может быть найден из условий нормирования. Наиболее распространен- 
ное условие нормирования -— равенство единице наибольшего отклонения от начального положе- 
ния. Форма поперечного колебания - это кривая изгиба оси балки в момент наибольшего откло- 
нения от нейтрального положения. При собственных колебаниях внешняя сила отсутствует, по- 
этому амплитуду колебаний определить невозможно, что и выражается в виде множителя с... 

При расчетах ряд (5) искусственно обрывают [1]. Ввиду этого, величины Р и &, являются 
приближенными, зависящими от количества членов ряда. С повышением использования количе- 
ства членов ряда точность расчетов растет. 
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Преобразованный вид формул для частот собственных колебаний дает следующее реше- 


ние 

2 _ Сы 2 Вы 

Ра а : +0, | ' (9) 

Ти Ти 

где 
К 
2 
Ти =: [тва 


0 


К 
в, = [5, © ав; 
0 
К 
Са — ЕТ, (5: ) : Аг, 
0 


И 
5 = [тьта» . 


В случае если ©, =0, то вращение не происходит, т.е. квадрат частоты собственных ко- 
лебаний невращающейся лопасти определяется по формуле 
2 Си 2 
РЕ Е: (10) 
Ти 


Для упрощения можно ввести коэффициент 


* 





Е, = вы : (11) 
Ти 
Соответственно формула (9) преобразуется к виду 
Ра = в № о, ы (12) 


Таким образом, собственные частоты колебаний лопасти можно определить, суммируя 
частоту собственных колебаний невращющейся лопасти и изменение частоты, вызванное увели- 
чением жесткости вследствие центробежной силы. Такой подход пригоден для определения соб- 
ственных частот колебаний лопасти. Коэффициенты т„, В», с„ вычисляются по формам собст- 


К 


венных колебаний невращающейся лопасти, также как Р 


квт 
обеспечить погрешность 1-2%. 

Выражение (12) показывает, что с ростом угловой скорости ветротурбины собственные 
частоты изгибных колебаний лопасти повышаются. Физически это объясняется действием расту- 
щей центробежной силы, стремящейся распрямить лопасть, т.е. увеличить ее жесткость на изгиб. 
Как известно, при повышении жесткости собственные частоты колебаний растут. 

Уравнение собственных колебаний невращающейся лопасти вытекает из уравнения (1) с 
учетом отсутствия центробежной силы 


. Согласно [1], данный метод может 


2 
(ЕТ „5”)"- р„т,ё=0. (13) 
При незначительном изменении удельной массы и жесткости по длине лопасти уравнение 
принимает вид 


Е” 8 =0, (14) 
2 
рт 
где а ==. 
ь + В 
Общее решение данного уравнения, выраженное через функции Крылова, имеет вид 
Е(+) = А5(аг) + ВТ(аг) + СИ (а') + РИ(аг), (15) 


где А, В, С, р - постоянные, определенные граничными условиями; 5, Т, (, Г - функции Крылова 
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5(о) = 0, 5(сВ а" + соза) 
Т(ог) = 0, 5(36 аи + эт о) 
И (ат) = 0,5(сВ ах — соз ог) 
Г (а) = 0,5(56 а" — зш аи") 
При г=0, 5(0) =Т, а Т(0) =У(0) =Ц(0) =0. Производные от одной из функция Крылова 
выражаются через другие: 


(16) 


о два Чссоъи = 0, 501 (58 си" — зто") = а, 


у. и г 
ЧТ 05 Вы = 0, 541 (св о" + соз аи) = 05, (17) 
Аг Аг Аг 
2 
ее Е. 
Аг Аг 


Таким образом, любая производная функции может быть найдена путем круговой замены. 

Для определения постоянных А, В, С, О используют четыре граничные условия [4]. Из 
(14) поочередно находятся производные и подставляются в соотношения при ’=0 и ’=^. Это дает 
систему уравнений относительно неизвестных А, В, С, О. Такая система имеет нетривиальное ре- 
шение. Тривиальное решение А=В=С=р=0 отвечает только состоянию покоя лопасти и не пред- 
ставляет интереса. Для нахождения нетривиального решения определитель, составленный из ко- 
эффициентов при неизвестных А, В, С, О, приравнивается к нулю. Так получают уравнение для 
нахождения с,, каждому корню которого соответствует свое нетривиальное решение исходной 


системы алгебраических уравнений. Подставляя последовательно каждое значение о, В систему 


алгебраических уравнений, определяются соотношения между коэффициентами 4, В, С, Р, т.е. 
соответствующую форму колебаний. Число собственных колебаний балок бесконечно, однако 
практический интерес представляет собой три-четыре значения аи,. 


Таким образом, частоты равны 





а (18) 


Приравняв определитель к нулю, получаем выражение для нахождения си, 








сВ(а.Ю) - соз(@К) +1=0. (19) 
Его решения равны 
1,571 
ок = } 
К 
_ 4,712 
ИЖ 9 
К 
20 
7,854 (20) 
о; = : 
К 
10,996 
и ==——. 
К 
Значение «а, может быть определено приближенно по формуле 
1-—0,5 
яд, 21 
: Е (21) 


Эта формула тем точнее, чем выше номер тона г. 
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Определение форм колебаний производится с учетом граничных условий винта 
5 (ог) (о,К) 
— Тею | 
Для построения форм колебаний целесообразно воспользоваться условиями нормирова- 
ния, в частности, равенства величины перемещения конца лопасти единице: 
К 
вю =ь [2 0)а"=1. (23) 
0 
Соответственно производные равны 
Е’(+) = Воь5 + СаТ + БаИ 
Е" (т) = Ва7У + Са? 5 + Ро?Т ь (24) 
Е" (0) = Во/У (0) + Са?5(0) + Ра?Т(0) = Са? —1=0 
Приведенная масса может быть выражена 


Е, (г) — С, И(а.) (22) 


74 у 
т, = [тёнаь = т| Енг =тА, (25) 
0 0 


В 
где А= [на 
0 


Численное интегрирование можно выполнить методом Симпсона: 
_ АС АС ЕС АС Е2С. РС 
о | 


Для вычисления ь, предварительно необходимо найти квадрат первой производной соб- 





А (26) 


ственной формы колебаний: 


* 
Б, — 


—-> 


( тт") (8' и а" = [= таг)’ та", (27) 


где С, =&, Дк - расстояние между расчетными сечениями лопасти. 


Резонансные параметры лопасти. На основании 
представленной методики был запрограммирован рас- 
чет колебаний лопастей малогабаритной ветроэнерге- 
тической установки, разрабатываемой ООО «Персо- 
нальные энергосистемы» (г. Ростов-на-Дону). Внешний 
вид установки представлен на рис.1. В расчете рас- 
сматривалась лопасть радиусом 2,25 м со стальным 
трубчатым лонжероном с жесткой заделкой в комле. 
Трубчатый лонжерон выполнял роль несущего элемен- 
та. К лонжерону крепился пенопластовый заполнитель 
с обшивкой. Согласно аэродинамическому расчету оп- 
тимальная частота вращения составила 380 об/мин. 
Расчет велся до частоты вращения 420 об/мин. Интег- 
рирование проводилось вдоль лопасти с шагом 0,1В. 

На рис.2 показаны расчетные формы первых 
трех тонов колебаний лопасти, приведенные по усло- 
вию нормирования к 1. 





Рис.1. Внешний вид малогабаритной 
ветроэнергетической установки 
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Рис.2. Формы собственных изгибных колебаний лопасти ветротурбины 


3-й тон 


ьная частота вращения 
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Рис.3. Резонансная диаграмма лопастей малогабаритной ветроэнергетической установки 


Как видно из резонансной диаграммы, представленной на рис.3, лопасти имеют хорошие 
резонансные характеристики, т.е. наступление резонанса возможно только на одной частоте, да- 
лекой от оптимальной. Для более наглядного представления значения по оси абсцисс представ- 


лены в относительных значениях, приведенных к оптимальной частоте вращения: 
=—, (28) 
Рот 


где р - оптимальная частота вращения ветроколеса. 
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На данном рисунке кривые тонов колебаний представляют собой графики собственных 
частот колебаний лопасти при увеличении центробежной силы. Наклонные прямые — графики 
частот вынужденных колебаний при различной частоте вращения ветроколеса при разных скоро- 
стях ветра. На рисунке видно, что первый тон собственных колебаний имеет низкие резонансные 
частоты, что характерно для всех ветротурбин, у которых лопасти имеют большой удельный вес и 
низкую жесткость. Пересечение первого тона с гармониками возбуждающей силы на частоте 
вращения менее 200 об/мин не представляет проблем, так как на этой частоте управление на- 
грузкой не осуществляется. 

Существует проблема с пересечением второго тона с пятой гармоникой на частоте 
297 об/мин. Для защиты от возможных колебаний в этой зоне необходимо добавить в алгоритмы 
управления требования или ослабления нагрузки, с целью выхода ветротурбины на частоты выше 
305 об/мин или увеличения, с целью понижения частот вращения ниже 290 об/мин. 

Проверить правильность выполнения расчета экспериментально можно, как правило, в 
одной точке — при отсутствии вращения. В этом случае на лопасти устанавливается виброгенера- 
тор малой массы и производится замер частоты появления резонанса по резкому росту амплиту- 
ды колебаний. Создать распределенное усилие, аналогичное нагрузкам на ветроколесе при вра- 
щении и сгенерировать колебания для нахождения резонансов представляется проблематичным. 
Выводы. Таким образом, представленный метод может быть использован для определения резо- 
нансных характеристик лопасти ветроустановок, имеющих малое изменение массы и жесткости 
вдоль лопасти. 

Расчет, выполненный для двухлопастной ветротурбины диаметром 4,5 м, показал, что не 
допустить появления резонанса на данном типе лопастей можно, использовав в алгоритмах 
управления нагрузкой ограничения рабочего диапазона частот 200-290 об/мин и 305-420 об/мин. 
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3. Мхеем В.А. Вазспуо* мепоуо{а па ргосНпоз(`. Ргосйпо$Е` |оразе] пези$ИНедо миа / В.А. 
Мхееу. -— М.: МАГ, 1973. - 205 $. — Тп Кизчап. 

4. ЭКортсем В.Г. ЗоБ$еппу`е коеБапуа юра$Н уепо!уока / В.Т. ЗКортсеу. -— Хаг` Ком: ХАТ, 
1983. - 200 $. — Тп Кизчап. 


ЗРЕСТЕУ1МС КЕЗОМАМСЕ ОТАСВАМ$ ОЕ СОМРАСТ АЕКОСЕМЕКАТОВ ВЕАОЕ$ 


А.$. СУВТМОМ, \.\У. РУОМТК 
(Роп За Тесйтса! Упмег®Ку) 


Тре гезопапсе спагацей$йс$ апг!уя5 {есптдие оЁ те сотрасЁ аегодепегаюг МааЕ 15 ргезетеа. Тйе аеюдп 
Юпт$ ое ПгзЕ [гее Маае ибгайоп юпеЕз апа а гезопапсе @адгат вауе Бееп БШ! Гог пе ипа игбте оЁ4.5 т 
т Фатег. 

Кеуигогаб: аегодепегайог, гезопапсе @адгат, ибгайоп тодеЕз. 
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